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Three-dimensional topological insulators are characterized by insulating bulk state and metallic surface 
state involving relativistic Dirac fermions which are responsible for exotic quantum phenomena and 
potential applications in spintronics and quantum computations. It is essential to understand how the Dirac 
fermions interact with other electrons, phonons and disorders. Here we report super-high resolution 
angle-resolved photoemission studies on the Dirac fermion dynamics in the prototypical Bi 2 (Te,Se) 3 
topological insulators. We have directly revealed signatures of the electron-phonon coupling and found that 
the electron- disorder interaction dominates the scattering process. The Dirac fermion dynamics in 
Bi 2 (Te 3 _ JC Se x ) topological insulators can be tuned by varying the composition, x, or by controlling the charge 
carriers. Our findings provide crucial information in understanding and engineering the electron dynamics 
of the Dirac fermions for fundamental studies and potential applications. 

Topological insulators represent a new state of matter that has an insulating bulk state but a conducting 
surface/edge state which hosts an unique spin texture and a linear dispersion relation compatible with 
massless Dirac fermions 1 " 3 . In addition to their potential applications in spintronics and quantum comput- 
ing 4 ' 5 , the topological insulators also provide a promising platform for realizing a number of novel and exotic 
quantum phenomena including quantum spin Hall effect 6 , fractional charge and quantized current 7 , dynamical 
axion field 8 , magnetic monopole 9 , Majorana fermions 1011 and quantum anomalous Hall effect 12 . Some of the 
phenomena have been observed experimentally in two-dimensional topological insulators, such as quantum spin 
Hall effect 13 and quantum Hall effect 14 in the HgTe quantum wells. However, in the case of three-dimensional 
topological insulators like Bi 2 (Te, Se) 3 family 15 " 18 , in spite of intensive efforts from transport experiments 19 " 23 , 
attempts to observe electronic transport signatures of the surface Dirac fermions have been hampered by a couple 
of obstacles. First, contrary to usual expectation that the bulk is insulating in topological insulators, the presence 
of defects in Bi 2 (Te, Se) 3 system makes the bulk conducting 24 and the transport properties are then dominated by 
the bulk state instead of the desired surface state 20 . Second, in several cases where the surface transport signatures 
are identified, more than one frequency or one type of charge carriers contribute to the quantum oscillation 21 " 23 , 
indicating the complexity of sample surface exposed to atmosphere or protection layers 25 . Third, the transport 
mobility of charge carriers in the Bi 2 (Te, Se) 3 topological insulators, fi' ~ 10 4 cm 2 V~ l s~ l (Refs. 19,23,26) is not 
sufficiently high compared to those from HgTe quantum wells(1.5 X 10 5 cm 2 W 1 ) 13,14 and graphene(2.3 X 
10 5 cm 2 V~ 1 s" 1 ) 27 " 29 . Although strict electron backscattering is forbidden in topological insulators due to time 
reversal symmetry invariance, it is important to uncover other scattering channels that limit the transport 
mobility of the two-dimensional Dirac fermions 26 . 

The charge transport properties are dictated by the electron dynamics in materials. In particular, the scattering 
of electrons by other entities such as other electrons, phonons and disorder, determines the lifetime of the 
quasiparticles, and thus the electron mobility and mean free path in the transport process 30 . Angle-resolved 
photoemission spectroscopy (ARPES) has become a powerful tool to directly probe such many-body effects 31 ' 32 . 
Under the sudden approximation, ARPES measures a single particle spectral function 



A(k,co) = 



\ImL(k,w) 
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with k and w representing electron momentum and energy, respect- 
ively. The real and imaginary parts of electron self- energy, ReY,(k, w) 
and Im£(/c, co), carry all the information about the band renormali- 
zation and quasiparticle lifetime in an interacting system 31,32 . 

There have been a few investigations on the dynamics of Dirac 
fermions in Bi 2 (Te, Se) 3 topological insulators using ARPES 33,34 or 
helium atom surface scattering probes 35 . However, the results 
obtained so far are not consistent. While the electron-phonon coup- 
ling strength (2) was reported to be 0.25 in Bi 2 Se 3 from temperature- 
dependent photoemission measurements 33 , it was claimed that the 
electron-phonon coupling is exceptionally weak in Bi 2 Se 3 in another 
measurement 34 . The inconsistency comes presumably from two 
issues related to the ARPES study on many-body effects in topo- 
logical insulators. First, the coexistence of bulk bands and surface 
bands near the Fermi level, due to the bulk doping by defects, com- 
plicates the identification of weak band renormalization effects in the 
surface bands. Second, the small momentum space and energy win- 
dow of the Dirac fermions and the weak many-body effects require 
high precision ARPES measurements. The momentum space occu- 
pied by the surface state, ~ ± 0.1 A" 1 for Bi 2 Se 3 (Fig. la) and ~ ± 
0.04 A" 1 for p-type Bi 2 Te 3 (Fig. 2a), is rather small 1718 ' 25 . To over- 
come these problems, it is ideal to first suppress the interference of 
the bulk bands in the measurements. Since the as-grown Bi 2 Se 3 is 
electron-doped (n-type) while Bi 2 Te 3 can be hole-doped (p-type) 25 , 
balancing the proportion of selenium and tellurium in Bi 2 (Se,Te) 3 
system may neutralize the bulk carrier doping and reduce the effect 
of the bulk bands on the surface state. Furthermore, vacuum ultra- 
violet (VUV) laser-based ARPES 36 can help suppress the bulk 
bands 25 , in addition to its super-high energy and momentum resolu- 
tions (Supplementary Fig. S2) that are desirable in studying many- 
body effects in topological insulators. 

Results 

Signatures of electron-phonon coupling in Bi 2 Se 3 . The Bi 2 Se 3 sam- 
ple shows sharp surface state band (Fig. la) that facilitates the 
investigation of the Dirac fermion dynamics in the material. The 
momentum distribution curves (MDCs) at different binding 



energies (Fig. lb) show well-defined sharp peaks that can be fitted 
by Lorentzian lineshape to obtain peak position and peak width 
(Full-Width-at-Half-Maximum, T MDC ) (Supplementary Fig. S3). 
The obtained MDC peak position as a function of the binding 
energy gives the measured dispersion relation (left panels of Fig. lc 
and Id). We assume a bare band as a linear line that connects the two 
points on the measured dispersion, one at the Fermi level and the 
other at 100 meV binding energy. This is reasonable within such a 
small energy window and particularly for the Dirac fermions that 
shows a linear dispersion over a large energy range. The effective real 
part of the electron self-energy (ReX) is obtained by subtracting the 
measured dispersion with the bare band, as shown in right panels of 
Fig. lc and Fig. Id (black empty circles). The imaginary part of the 
electron self-energy (ImX) is obtained from the fitted MDC width: 
ImE = r MDC X v 0 /2 with v 0 being the Fermi velocity of the bare band 
(right panels of Fig. lc and Fig. Id, blue empty squares) 32 . 

The extraction of the electron self-energy provides full informa- 
tion for understanding the electron dynamics in a material where 
the electron scattering usually involves three main processes: 
electron- electron interaction, electron-phonon coupling, and elec- 
tron-disorder interaction 37,38 . The electron-disorder interaction is 
characterized by the imaginary part of the electron self- energy at 
zero energy (Fermi level) since the contribution from the electron- 
electron interaction and electron-phonon coupling approaches zero 
at the zero energy at very low temperature. The electron-phonon 
interaction is characterized by the band renormalization (usually a 
kink in dispersion, corresponding to a peak in the real part of the 
electron self-energy) and a drop in the imaginary part of the electron 
self-energy. As exemplified in Fig. 1 for Bi 2 Se 3 , there are clear indica- 
tions of electron-phonon coupling revealed, such as the kink in dis- 
persions (left panels of Fig. lc and Fig. Id), peak in the real part of the 
electron self- energy and a drop in the imaginary part of the electron 
self- energy (right panels of Fig. lc and Fig. Id). The characteristic 
energy scale of the phonons involved (co 0 ) is determined by the peak 
position of the real part of the electron self- energy; in the case of 
Bi 2 Se 3 , it is — 18 meV as marked by arrows in Fig. lc and Id. The 
electron-phonon coupling strength, 2, can be determined by the 




Figure 1 | Observation of electron-phonon coupling in Bi 2 Se 3 . (a) Surface state dispersions along T — K(\eft panel) and T — M(right panel) directions. 
The inset shows the Fermi surface of Bi 2 Se 3 . (b) Momentum distribution curves (MDCs) of the surface state dispersion along T — K. Red line and blue line 
indicate the MDCs at the Dirac point (E D ) and the Fermi level (Ep), respectively, (c) Dispersion and effective electron self-energy of surface state 
band along T — ^direction. The left panel shows the measured dispersion (red circles) obtained from MDC fitting. The blue dotted line represents an 
empirical bare band that is a straight line connecting two points in the measured dispersion, one at E F and the other at 100 meV binding energy. The right 
panel shows the effective real part of the electron self-energy, ReL{k, co) (black circles), and imaginary part of the electron self-energy, ImL(k y co) (blue 
squares). The red dashed line is a linear fit to the low-energy region of ReL(k, co). Blue arrows indicate the characteristic energy scale of the electron- 
phonon coupling. (d)Dispersion and effective electron self-energy of surface state band along T — Mdirection. 
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Figure 2 | Extremely weak electron-phonon coupling in p-type Bi 2 Te 3 . (a) Surface state dispersions along T — K(\eft panel) and T — M(right panel) 
directions. The inset shows the Fermi surface of p-type Bi 2 Te 3 . (b) Momentum distribution curves (MDCs) of the surface state dispersion along T — K. 
Red line and blue line indicate the MDCs at the Dirac point (E D ) and the Fermi level (E P ), respectively, (c) Dispersion and effective electron self-energy of 
surface state band along T — ^direction. The left panel shows the measured dispersion (red circles) obtained from MDC fitting. The blue dotted line 
represents an empirical bare band that is a straight line connecting two points in the measured dispersion, one at E F and the other at 40 meV 
binding energy. The right panel shows the effective real part of the electron self-energy, ReL(k y w) (black circles), and imaginary part of the electron 
self-energy, ImL(k, co) (blue squares). (d)Dispersion and effective electron self-energy of surface state band along T — Mdirection. 



slope of the line near the zero energy in the real part of the electron 
self-energy. As shown by red dashed line in the right panel of Fig. lc, 
it is determined to be 2—0.17 for Bi 2 Se 3 . The electron -electron inter- 
action manifests itself more clearly in the imaginary part of the 
electron self- energy. In a typical Fermi liquid case, its contribution 
is proportional to the (E-E F ) 2 : lm ee = /?(E-E F ) 2 . Since the electron- 
phonon contribution is constant above the characteristic phonon 
scale, the deviation from the constant above the phonon scale gives 
a measure of the electron-electron interaction involved. For Bi 2 Se 3 
and the Bi 2 (Se,Te) 3 series, we find that the electron -electron contri- 
bution is weak (/?— 0.014 for Bi 2 Se 3 ), particularly when compared 
with the dominant electron-disorder interaction and the electron- 
phonon coupling. Therefore, in the following, we will mainly focus 
on the electron-phonon coupling and electron- disorder interaction. 

As shown in Fig. 1, our measurements clearly indicate that there is 
an obvious electron-phonon coupling (2—0.17) in the electron 
dynamics of the Dirac fermions in Bi 2 Se 3 . It is nearly isotropic 
because similar behaviors are observed along the f — K(Fig. lc) 
and r — M(Fig. Id) momentum cuts. The characteristic energy scale, 
— 18 meV, is consistent with the highest-frequency phonon mode 
(Aj 2 g ) at — 20 meV in Bi 2 Se 3 from both theoretical calculations 39 and 
Raman measurements 40 (Supplementary Fig. S7). This is in strong 
contrast to a recent report where exceptionally weak electron-pho- 
non coupling was claimed and no dispersion kink was resolved in 
Bi 2 Se 3 34 . We believe that the clear manifestation of electron-phonon 
coupling in our measurements is most likely due to much improved 
instrumental resolution we have achieved. 

Extremely weak electron-phonon coupling in p-type Bi 2 Te 3 . The 

electron-phonon coupling is extremely weak in the p-type Bi 2 Te 3 
sample, as seen in Fig. 2. This sample was prepared under similar 
condition as the above Bi 2 Se 3 sample and it is p-type because the 
Fermi level intersects with the bulk valence band 25 . Even measured 
with such a super-high instrumental resolution (Supplementary Fig. 
S2), we do not resolve any clear signature of electron-phonon 
coupling. The measured dispersion is basically a straight line, 
corresponding to a real part of electron self- energy that is nearly 
zero (black empty circles in right panels of Fig. 2c and Fig. 2d). 



There is no clear drop in the imaginary part of the electron self- 
energy (blue empty squares in right panels of Fig. 2c and Fig. 2d). 
All these observations point to an extremely weak electron-phonon 
coupling in the p-type Bi 2 Te 3 sample. In the mean time, the p-type 
Bi 2 Te 3 also exhibits a rather weak electron- disorder scattering when 
compared with Bi 2 Se 3 . As seen from the right panel of Fig. 2c, the 
imaginary part of the electron self- energy at zero energy is only 
5.6 meV along the T — ^direction. This is significantly lower than 
37 meV along the same direction in Bi 2 Se 3 topological insulator 
(right panel of Fig. lc). 

Tunability of Dirac fermion dynamics in the Bi 2 (Te 3 _ x Se x ) series. 
The dramatically different behaviors between Bi 2 Se 3 (Fig. 1) and 
p-type Bi 2 Te 3 (Fig. 2) motivated us to further investigate the com- 
position dependence of the Dirac fermion dynamics in the 
Bi 2 (Te 3 _ x Se x ) series. We observe a systematic evolution of band 
structure, electron-phonon coupling, and electron-disorder 
scattering, with the change of composition x in Bi 2 (Te 3 _ x Se x ) 
(Fig. 3). Note that all these samples are prepared under similar 
conditions. First, the location of the Dirac point (E D ) moves 
monotonically to higher binding energy with increasing x. For 0 < 
x < 0.05, the samples are p-type because the Fermi level intersects 
with the bulk valence band. For 0.15 < x < 1.5, the Fermi level 
intersects with the surface state bands only. For 1.5 < x < 3, the 
samples become n-type because the Fermi level intersects with the 
bulk conduction band. Note that the Dirac point exhibits a big jump 
when x changes from 0.05 to 0.15. Second, there is an overall increase 
of electron-disorder interaction with increasing x in Bi 2 (Te 3 _ x Se x ). 
This can be seen from the increase of the imaginary part of the 
electron self- energy and hence its value at zero energy with 
increasing x. This indicates that the Bi 2 (Te 3 _ x Se x ) system gets 
more disordered with increasing x. Third, there is also an overall 
increase in the electron-phonon coupling strength with increasing 
x in Bi 2 (Te 3 _ x Se x ). Although no indication of electron-phonon 
coupling is observed in p-type Bi 2 Te 3 , a slight substitution of Te by 
Se, with x = 0.05, already renders a clear observation of electron- 
phonon coupling in Bi 2 (Te 2 .9 5 Se 0 .o5)> as seen in the appearance of a 
peak in the real part of the electron self- energy (black circles in 
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Figure 3 | Evolution of surface state band structure and electron self-energy with the composition, x, in Bi 2 (Te 3 _ x SeJ system, (a) The surface state 
band structure of Bi 2 Te 3 -^>e x with various compositions, x, measured along the F — ^direction, (b) The surface state band structure of n-type 
Bi 2 Te 3 along the F — ^direction, (c) Corresponding effective real part of the electron self-energy (black circles) and imaginary part of the electron 
self-energy (blue squares) for Bi 2 (Te 3 _ x SeJ with different compositions, (d) Corresponding effective real part of the electron self-energy (black circles) 
and imaginary part of the electron self-energy (blue squares) for n-type Bi 2 Te 3 . 



Fig. 3c2). The electron -phonon coupling strength increases from 
nearly zero for p-type Bi 2 Te 3 to 2—0.17 for Bi 2 Se 3 . 

It is interesting that for the same Bi 2 Te 3 nominal composition, 
when prepared under different conditions, samples can become 
either p-type or n-type that exhibit dramatic difference in their 
Dirac fermion dynamics. As shown in Fig. 3b, the Dirac point of 
the n-type Bi 2 Te 3 lies nearly 310 meV below the Fermi level which 
differs significantly from the 85 meV location in the p-type Bi 2 Te 3 
(Fig. 3al). Particularly, the n-type Bi 2 Te 3 exhibits a pronounced 
electron-phonon coupling, as seen from the obvious peak in the real 
part of the electron self- energy (black circles in Fig. 3d) and a drop in 
the imaginary part of the electron self- energy (blue squares in 
Fig. 3d). The resultant electron-phonon coupling constant is 
A — 0.19, which is even slightly larger than that in Bi 2 Se 3 (2—0.17). 
The electron-disorder scattering also increases in the n-type Bi 2 Te 3 
compared to the p-type one, as the value of the imaginary part of the 
electron self energy at zero energy becomes much larger (Fig. 3d). 
However, compared with Bi 2 Se 3 , the disorder scattering in n-type 
Bi 2 Te 3 remains much weaker(Fig. 3cll). The first principle calcula- 
tion, which expects a rather weak electron-phonon coupling in 
Bi 2 Te 3 (2— 0.05) 41 , is hard to explain the strong electron-phonon 
coupling in n-type Bi 2 Te 3 and particularly the disparity of elec- 
tron-phonon coupling between the n-type and p-type Bi 2 Te 3 . 

The different electron dynamics between Bi 2 Se 3 and Bi 2 Te 3 , and 
between n-type Bi 2 Te 3 and p-type Bi 2 Te 3 , is intriguing. We note that 
the native majority defect types are totally different in these samples, 
i.e., Se vacancies in Bi 2 Se 3 , Te-Bi antisites in the n-type Bi 2 Te 3 , and Te 
vacancies in the p-type Bi 2 Te 3 . The difference of the electron- dis- 
order interaction can be due to such different types of disorders and 
their densities in the samples. In addition, as the main difference 
between the n-type and p-type Bi 2 Te 3 lies mainly in the relative 



position between the Fermi level and the bulk bands, it is reasonable 
to deduce that the bulk bands may play a role in the electron 
dynamics of the Dirac fermions. This is consistent with the obser- 
vation that, for n-type Bi 2 Te 3 (x = 0), although its composition is 
rather different from Bi 2 (Te 2 85 Se 0 15 ) (x = 0.15), their Dirac fermion 
dynamics is quite similar in both the electron-phonon coupling and 
electron-disorder interaction (Figs. 3 and 4). The distinction of the 
electron-phonon coupling between the n-type and p-type Bi 2 Te 3 
indicates that the control of the charge carrier doping can be an 
alternative way to tune the Dirac fermion dynamics in topological 
insulators. 

One may wonder whether the origin of the observed inelastic 
scattering in Bi 2 (Se,Te) 3 -based topological insulators is related to 
the materials themselves (defect, atomic arrangement and etc.) or 
the intrinsic topological properties associated with the time reversal 
symmetry (TRS) invariance. The TRS-protection may reduce the 
scattering strength because strict backscattering is not allowed, but 
it is difficult to quantitatively analyze this weakening effect. Note that 
the entire Bi 2 (Te,Se) 3 series are topological insulators. Therefore, we 
believe the systematic evolution of their electron dynamics is mainly 
a manifestation of the usual many-body effects, other than a mani- 
festation of the intrinsic TRS-related topological properties. 

Fig. 4 gives a quantitative summary of the electron-phonon coup- 
ling and electron-disorder scattering in the Bi 2 (Te 3 _ x Se x ) system. 
The relative position between the Dirac point and the bulk band is 
depicted in Fig. 4a for a couple of typical compositions and both p- 
and n-type Bi 2 Te 3 . With increasing x in Bi 2 (Te 3 _ x Se x ), there is an 
abrupt jump in the electron-phonon coupling strength from x = 0 to 
x = 0.15, accompanied by a transition from p-type to n-type. Further 
increase of x leads to an overall slight increase of electron-phonon 
coupling (Fig. 4c). Meanwhile, the characteristic energy scale, co 0 , 



SCIENTIFIC REPORTS | 3 : 241 1 | DOI: 1 0.1 038/srep0241 1 



4 




Bi 2 Te 3 (p) 



Bi 2 Te 3 (n) 



Figure 4 | Variation of electron-phonon coupling and electron-disorder scattering with composition, x, in Bi 2 (Te 3 _ x SeJ topological insulators. On 

the far right side, data on n-type Bi 2 Te 3 are also included, (a) Schematic band structure at several typical compositions to show the relative position 
between the topological surface states (color dashed lines) and the bulk states (green solid lines). Black dashed lines indicate the Fermi level, (b) Variation 
of the electron-disorder scattering (red circles) with composition, x, in Bi 2 (Te 3 _ x Se x ) topological insulators. It is obtained by taking the value of the 
imaginary part of the electron self-energy at zero energy (E P ) as shown in Fig. 3c. Also shown is the electron mean free path (blue circles) calculated from 
the measured MDC width at the Fermi level, (c) Variation of the electron-phonon coupling strength X (red squares) and characteristic energy scale co 0 
with the composition, x, in Bi 2 (Te 3 _ J£ Se x ) topological insulators and in n-type Bi 2 Te 3 . The dashed lines in (b) and (c) are guide to the eyes. 



also increases with x, which is consistent with the increase of the 
phonon frequency from Bi 2 Te 3 to Bi 2 Se 3 (Supplementary Fig. S7). 
This indicates that the (Te,Se) vibrations play a major role in dictat- 
ing the electron-phonon coupling strength in the Bi 2 (Te 3 _ x Se x ) 
system. This is also consistent with the observation that the char- 
acteristic energy scale is identical in n-type and p-type Bi 2 Te 3 because 
their phonon frequencies are quite similar (Supplementary Fig. S7) 
although they have different electron-phonon coupling strength. It is 
clear that electron -disorder scattering is dominant in the scattering 
process. As seen from Fig. 4b, the imaginary part of the electron self- 
energy at zero energy Im£(E F ), which gives a good measure of the 
electron scattering rate, shows an overall increase from the p-type 
Bi 2 Te 3 for x = 0 to Bi 2 Se 3 for x = 3 in the Bi 2 (Te 3 _ x Se x ) system (red 
circles in Fig. 4b). The quasi-particle mean free path, which can be 
determined from the MDC width r MDC ,€' = l/r MDC , first decreases 
abruptly from the p-type Bi 2 Te 3 (x = 0) to x = 0.15, followed by an 
overall slight decrease with x in Bi 2 (Te 3 _ x Te x ) (blue circles in 
Fig. 3b). The electron mobility of the Dirac fermions (//) is directly 
related to the quasi-particle mean free path, \i' = e£'/(hk F ), with k F 
being the Fermi momentum. It is then estimated that the electron 
mobility in p-type Bi 2 Te 3 is near 1000 cm 2 V^s -1 (Supplementary 
Fig. S4). By comparison, the electron mobility in Bi 2 Se 3 and n-type 
Bi 2 Te 3 is ~ 70 and —100 cm 2 V~ 1 s~ 1 , respectively, which are signifi- 
cantly lower than that of p-type Bi 2 Te 3 . We note that, although the 
electron mobility obtained from ARPES is not strictly the same as the 
electron mobility in the transport properties, they have good corres- 



pondence, especially for materials with circular Fermi surface and 
nearly isotropic electron scattering. 

Discussion 

The present study will provide important information for engineer- 
ing the Dirac fermion dynamics in Bi 2 (Te 3 _ x Se x ) topological insu- 
lators. It clearly indicates the presence of electron-phonon coupling 
and dominant role of the electron-disorder scattering in dictating the 
Dirac fermion dynamics in the system. The Dirac fermion dynamics 
can be tuned by varying either the composition or the charge carrier 
concentration. Specifically for the Bi 2 (Te 3 _ x Se x ) system, our results 
clearly show that Bi 2 Te 3 system is superior over Bi 2 Se 3 in achieving 
higher electron mobility both due to its weak electron-phonon coup- 
ling and its weak electron -defect scattering. To further enhance the 
electron mobility in searching for novel quantum phenomena, it is 
necessary to choose a system with weak electron-phonon coupling 
and reduced disorder scattering. 

Methods 

High quality single crystals of Bi 2 (Te 3 _ x SeJ (0 < x < 3) were grown by the self-flux 
method. The mixed materials were heated to 1000°C, held for 12 hours and then 
slowly cooled down to 500°C over 100 hours before cooling to room temperature. 
Single crystals of nearly one centimeter in size were obtained by cleaving. The angle- 
resolved photoemission measurements were carried out on our vacuum ultra-violet 
(VUV) laser-based angle- resolved photoemission system, with an overall energy 
resolution of ~ 1 meV. The momentum resolution is ~ 0.004 A -1 for 30° angular 
mode and ~ 0.002 A -1 for 14° angular mode of the R4000 electron energy analyzer. 
The samples were all cleaved and measured in vacuum with a base pressure better 
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than 5X10 11 Torr and a temperature T ~ 20 K. For more details about the 
experimental methods, see supplementary information. 
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